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ABSTRACT 

The origin of the long secondary periods (LSPs) in red variables remains a mystery up to 
now, although there exist many models. The light curves of some LSPs stars mimic an eclipsing 
binary with a pulsating red giant component. To test this hypothesis, the observational data 
of two LSP variable red giants, 77.7795.29 and 77.8031.42, discovered by the MACHO project 
from the LMC, are collected and analyzed. The probable eclipsing features of the light curves 
are simulated by the Wilson-Devinney (W-D) method. The simulation yields a contact and a 
semidetached geometry for the two systems, respectively. In addition, the pulsation constant of 
the main pulsating component in each binary system is derived. By combining the results of the 
binary model and the pulsation component, we investigate the feasibility of the pulsating binary 
model. It is found that the radial velocity curve expected from the binary model has a much 
larger amplitude than the observed one and a period double the observed one. Furthermore, the 
masses of the components based on the density derived from the binary orbit solution are too 
low to be compatible with both the evolutionary stage and the high luminosity. Although the 
pulsation mode identified by the pulsation constant which is dependent on the density from the 
binary-model is consistent with the first or second overtone radial pulsation, we conclude that 
the pulsating binary model is a defective model for the LSP. 



Subject headings: binaries: close - galaxies: indiv 

1. INTRODUCTION 

Among the variable red giant stars, one sub- 
type exhibits long secondary periods (LSPs). The 
light curves of these stars exhibit not only a short 
primary period but also a long secondary period, 
which is approximately nine times longer than the 
short one. This phenomenon has been known 



idual (LMC) - stars: late-type - stars: oscillations 

for se veral decades ( Pavne-Gaposhkiiil 1954 ; Houk 



19631 ) . So me sample s of th ese LSPs variables are 



shown in Kiss et al 



(|l999l ). An interest in the 
st ars with the LSPs h as been renewed by the study 
of I Wood et all (Il999l) . This paper shows that in 
the LMC, ~25 % of all variable asymptotic gi- 
ant branch (AGB) stars show LSP roughly nine 
times longer than the short primary period which 
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is typically ~ 30-200 days. Meanwhile, a study 
of the bright local pulsating red giants indicates 
that at least one-th ird of these stars exhibit LSPs 
(jPercv et alj 120041) . Soszynski also gives ~30% 



Wood et al.ll2004 iNichoUs et al.ll2009l) . Some ob- 
servational arguments supporting the binary hy- 



of pulsating red giants in the LMC with LSPs 
(jSoszynski et al.l 120071. In the period- luminosity 



(P-L) diagram, it is interesting to see that the LSP 
variables follow a distinct sequence (sequence D), 
which is roughly parallel to the radial pulsation 
sequences A, B, and C for variable red giants. 

The LSPs present in variable red giants have 
attracted a lot of attention since their discovery, 
but their origin still remains mysterious. Since the 
LSPs are several times longer than the fundamen- 
tal radial periods, they could not b e caused by 



normal ra dial pulsations. Mo reover, iWood et al 



(|l999l ) and I Wood eraP (|2004 ) note that the LSPs 
can not be explained by mode for the oscilla- 
tory mode is evanescent in convective region 
and it is unlikely to be observable in a red gi- 
ant because of its convective envelope. The g~ 
mode is dynamically unstable in the convection 
envelope, and unable to lead to any oscillation. 
Regarding the nonradial p modes, their periods 
are rather shorter than those of the fundamental 
radial modes, so they can not explain the LSPs 
either. Alternatively, if the LSPs are caused by 
some strange modes, they wo uld be extr emely 
damped and should not be seen (|Woodll2000l) . One 
more possible explanation is the rotating spheroid 
model, which can explain the shape of the veloc- 
ity curve, but there is no reason for the rotating 
period to bring about the observed P-L relation. 

It seems that variable red giants with long sec- 
ondary periods (hereinafter referred to as LSPVs) 
can not be easily interpreted as pulsating red gi- 
ants. In this situation, a hypothesis of bina- 
rity arose. It has been suggested that the se- 
quence D stars could be components of close bi- 
nary systems, and the LSPs could be interpreted 
as the light variations caused by ellipsoi d al bi- 



nary motions or eclipses. Soszvhski et al.l (|2004 , 



20071 ) find that the sequence D stars overlap with 
and have a direct continuation of the sequence 
E stars that are mostly confirmed to be bina- 
ries. Radial velocity variations on the time-scale 
of LSPs have also been measured in a number of 
sequence D stars, such radial velocity variations 
being a requirement if the sequence D stars are 
binaries (|Hinkle et al.ll2002l : lOlivier fc WooJ2003t 



pothe sis have already been reviewed bylSos zvhski 
(I2OO7I) 



After examining the light curves of LSPVs col- 
lected from the MACHO database, we note that 
many of them show some nearly regular and sta- 
ble, eclipse-like light variations at the long sec- 
ondary periods with a large amplitude in com- 
parison with that of the primary pulsation. It 
seems that the LSPs can be easily interpreted as 
eclipses by an orbiting component. If this can 
be shown to be the case, it would be direct ev- 
idence to support the binary hypothesis for the 
sequence D stars. To test this idea, we propose 
a model - pulsating binary, i.e. a binary system 
with a pulsating component. We begin with some 
very probable eclipsing LSPs candidates whose 
light curves look like eclipsing binary with pri- 
mary and secondary minima. The light curves of 
two LSP stars are collected and analyzed by us- 
ing the Wilson-Devinney ( W-D, hereafter) code 
and the power-spectrum method. Afterwards, the 
proposed eclipsing and pulsating nature as well as 
their evolutionary properties are discussed. 

2. THE DATA 

Since the LSPs are sometimes as long as 1500 
days, the candidate selection is based on the long- 
term MACHO survey The MACHO project is 
a microlensing experiment that monitored numer- 
ous stars in the LMC, SMC and Galactic Bulge 
over a time of ^ 3000 days. As the by-products of 
this survey, a large number of variables were dis- 
covered and monitored, including many LSPVs. 
From this database, two LSPVs in the LMC are 
selected as our program objects. Their F.T.S 
(field. title. sequence) numbers are 77.7795.29 and 
77.8031.42, respectively. They are chosen as the 
working sample for the following reasons: (1) The 
time base of observations for both of the stars is 
long enough to cover at least a couple of LSPs. 
(2) Their light variations are both eclipse-like and 
symmetric over the long secondary period, and 
have distinct short primary variation. (3) These 
two stars are quite different from each other not 
only in their LSPs but also in the shapes of the 
light curves, they are expected to be different 
types of binaries. 
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The MACHO data are taken in non-standard 
red and blue bandpasses. The adopted trans- 
formation from the MACHO instrumental pho- 
to metry to Kron-Cou sins R and V are given 



Alcock et al 



(I1999I) . However, it does not 
work well for our candidates in practice and 
brings about visible dispersion for the photometric 
points. Analyzing the data sets, we find the reason 
mainly being that the observed c olor i ndex term 
"(VM,t - RM,tT in I Alcock etldl (|l999( ) brings in 
large dispersion as the raw red and blue magni- 
tudes transfer their errors to each other. This 
term requires both measurements to be good at 
each time when adopting the transformation. In 
addition, the transformation has a term involving 
airmass. When the value of the airmass is greater 
than 2.0, the correction for the atmosphere ex- 
tinction is rather uncertain. We adopt another 
transformation to a void the above problems as 
Struble et all (|2006l) : 



R = Rm.i + 23.90 + 0.1825(14/,* - Rm.i) (la) 
V = Vma + 24.22 - 0.1804(14/.* - Rmx), (lb) 

where the term "Va/,* — RM,t^ is the mean color 
index for a red giant instead of an individual mea- 
surement. Since color indexes of the red giants 
have a small scattering of about 0.2 magnitude, 
the adoption of the mean multiplied by a factor of 
0.18 brings about a dispersion of only 0.04 mag, 
much smaller than the photometric uncertainty. 
The light curves of standardized R and V magni- 
tudes for our two candidates are shown in Fig. 1. 

OGLE-II ()Udalski et al.ll997l:ISzvmahskil2005[) 
also provides the light variation data for red giants 
in the I band, and the time span of the light curve 
overlaps the last half of the MACHO data, with 
~ 300 days extension which is not long for the 
LSP. Besides and more importantly, the data are 
much sparser than the MACHO in the overlapping 
time range and it is of little help in analyzing the 
short primary pulsation. Thus we don't include 
the OGLE-II data for following analysis. 

3. THE BINARY MODEL 

Assuming that the phenomenon of the LSPs 
is caused by an orbiting component, we can em- 
ploy some binary simulation method to analyze 
the LSPs light curve. The W-D method, which 
consists of two main FORTRAN programs - LC 



and DC (short for Light Curve and Differential 
Correction, respectively), is a simulation method 
producing the photometric solution for the binary 
system. 

The 2007 version of the W-D code is used to 
analyze the eclipse-like light curves in R and V 
bandpasses of our two candidates. Nonlinear limb- 
darkening via a log arithmic form along wi t h many 
other featu r es (iWilson &: Devinnevl llQTlt IWilson 
19791 11990I: iKallrath et al.lll998l ) are used in the 
code. Considering the likely close distance be- 
tween the components, the effect of reflection is 
taken into account. 

In computing the photometric solution, the im- 
portant parameters adopted in the DC program 
are as follows. The orbital period for adjustment 
is computed by the PDM (Phase Disp ersion Mini- 
mization) method ( Stellingwer j|l978[ ) and the Pe- 
riod04 software for double-check at the beginning. 
The value is about 970 days for 77.7795.29 and 
1700 days for 77.8031.42. Note that these or- 
bital periods are twice the normally adopted value 
for the LSP, the latter corresponding to one cy- 
cle of long period variation apparent in the light 
curve. The temperatures of the two primary stars 
are bo th set to 3311 K according to IFox &: Wood 
(Il982h for long period stars. The initial bolometric 
(X1,X2,Y1,Y2) and monochromatic (xl,yl,x2,y2) 
limb-darke ning coefficients of the components are 
taken from lVan Hamm^ (1993). The gravity dark- 
ening exponents are set to 0.32 for both the pri- 
mary and the secondary comp o nent for convective 
envelopes according to iLucvl (Il967h . The bolo- 
metric al bedos are taken as Al=0.4 and A2=0.5 
following iRucinskil (|l969t) . 

The preliminarily estimated parameters are put 
in the first iteration of the DC code and after sev- 
eral iterations, a converged solution is reached. 
The photometric solutions are shown in Table 1, 
and the synthesized R band light curves computed 
from the LC code are shown in Fig. 1 together 
with the observed curves of the two candidates. 
The photometric solutions in the table provide in- 
sight into the orbital dynamics of the binary sys- 
tem. The orbit period (Porb) adjusted by the DC 
code is consistent with the value from the PDM 
and Period04 method. The large angles of inclina- 
tion (i), and argument of periastron (w) describe 
the orientation of the binary orbit, while the value 
of e reveals a nearly circular orbit. Moreover, we 
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find that the secondary component is a low mass 
M-type star, according to the mass ratio and effec- 
tive temperature. All the binary modes, including 
detached mode, semidetached mode, and contact 
mode in the code are tried in the solution-seeking 
procedure, and it is found that only one mode 
can reach a converged solution for each candidate. 
For 77.7795.29, a contact configuration converges, 
and for 77.8031.42, a semidetached system with 
the primary component filling its Roche lobe con- 
verges. This implies that both of them have a mass 
exchange with their components due to Roche lobe 
overflow. No converged solutions were found for 
systems showing purely ellipsoidal variations i.e. 
systems with an invisible secondary companion. 

Using the photometric solutions, the theoret- 
ical radial velocities can be also simulated by 
the LC code. We try this in order to compare 
them with some observational facts. The bot- 
tom panel of Fig. 1 shows the synthesized ve- 
locity curves for the primary components. The 
full amplitude of the velocity variation is about 
11.0 km s-i for 77.7795.29, and 9.0 km s'^ for 
77.8031.42, where the value of the binary system 
mass is taken from later Section 5 (but see discus- 
sion there for the problem of mass). Yet all the 
studi e s of radial velocity for L SPVs (iHinkle et al 



20021: lOlivier fc WoodI 120031: IWood et all 120041: 

NichoUs et al.l l2009h show typical amplitudes of 
3-4 km s~^ and no objects with full velocity am- 
plitude greater than 7 km s~^ are found. More- 
over, the synthesized velocity curves of the pri- 
mary components have a period which is twice 
of the light variation period. In contrast, all the 
observed radial velocities for the LSPVs have the 
same period as that of the light variation. The 
synthesized radial velocities for the second stars 
are also obtained from the LC code, and we find 
that the full amplitude is about 18.5 km s~^ for 
77.7795.29 and 18.3 km s^^ for 77.8031.42. The 
large velocity separations between the primary 
and secondary stars would have been readily seen, 
while in the existing spectral observations there 
is no indication of this. Therefore the synthe- 
sized velocity curve based on the binary model is 
inconsistent with the observation. 

In addition, we should note that, among all the 
parameters derived from the DC code, some are 
highly reliable, including the orbital period, mass 
ratio and temperature difference if we admit the 



binary hypothesis, and some are relatively uncer- 
tain, such as the effective temperatures of the two 
components. Several parameters derived from the 
LC code, such as the mass and radius, are par- 
ticularly uncertain for there are no radial veloci- 
ties available for comparison, so they are not listed 
in Table 1. The synthesized radial velocities pro- 
duced by the LC code are also very uncertain. 

4. THE INTRINSIC PULSATION 

For an eclipsing binary with a pulsating (pri- 
mary) component, the observed light curve could 



be a pproximately interpreted as (jZhang et al 
2009h : 

^obs = hfpui + h, (2) 

where h and I2 are the calculated brightness con- 
tributed by the primary and secondary compo- 
nents respectively, and /pui denotes the pulsating 
variation of li. With the derived binary photo- 
metric solution, we could compute the brightness 
of the binary components separately at each ob- 
servation epoch by using the LC code in the W- 
D code package, from which the "pure" pulsation 
light variation of the primary star could then be 
extracted. However, we do not use this approach, 
because the intrinsic pulsating amplitude of the 
primary component is so large that the theoreti- 
cally synthesized light curve computed by the LC 
code would have up-down fiuctuation and the ex- 
tracted light curve would not be exactly equivalent 
to the true intrinsic pulsation variations. Instead, 
we make use of the original data sets to analyze 
the primary pulsation. 

Period04, a commonly used technique for time 
series analysis, which utilizes Fourier transforms 
as well as multiple-least-squares algorithms, is ap- 
plied to the power spectrum analysis. In doing 
that, we use the R band time series data which 
are brighter and with higher photometric accu- 
racy than in the V band. We select only those 
peaks in the power spectrum with signal to noise 
ratio (S/N) larger than 4.0 for further discussion. 
Fig. 2 and Fig. 3 show the first four-frequency so- 
lution of the data sets in the R band before and 
after subtraction of the most prominent remaining 
frequency (prewhitening) . At the top of the two 
figures the spectral window based on the epoches 
of available observations is displayed. We notice 
that the alias patterns, including the 1 c/d daily 
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alias, are quite low in power. The next four pat- 
terns in both figures are the step-by-step power 
spectrum derived from the pulsation data. The 
fitting is stopped when the light curve is well fit- 
ted. For 77.7795.29, nine frequencies are obtained, 
by using most parts of the MACHO data while the 
last part of the MACHO data is dropped because 
of a large gap. Among all the nine frequencies 
obtained, we find that the main power spectrum 
are dominated at /i = 0.0020 c/d and /a = 0.0100 
c/d. /i is two times the assumed orbital frequency 
(/o = 1/970 = 0.00103 c/d), it is not accepted as 
the real pulsating frequency, so only /2 remains as 
the intrinsic pulsating frequency. The other de- 
rived frequencies are found to be related to /o and 
/2 as follows: /i = 2/o, /s = /2 + 6/o, h = /o, 
/5 = 4/o, /e = /2 + 1/365, /t = /2 + 7/0, /s = 3/o, 
/g = 5/0, all of them are aliases. For the other 
star 77.8031.42, we obtain /2 = 0.0112 c/d, the 
real pulsating frequency, and the others have the 
relations of: /i = 2/o, h = 4/o, U = f2 + 0.0004, 
/5 = 3/6, /e - 2/0 + fs, h = 4/6, /s = 1/365, 
/g = /2 + 10/0, /lo - /2 - 0.0004, /n = /6 + /g, 
/12 = 10/0, /i3 12/0, /i4 = /o, where /o is the 
assumed orbital frequency. Here, 0.0004 c/d, may 
be the rotation frequency of the LSPVs, in the the- 
oretical range of 240 0-10000 days for the rotation 
perio d of red giants (jKiss et al. I I2OOOI: IWood" et al 



20041 ). Based on the above analysis, only one short 
primary pulsation is found for our two candidates 
respectively an d the result is consistent with previ- 
ous result from Wood et"al] ( 19991) for LSPs stars. 

The eclipsing light curve synthesis has provided 
us some important parameters of the binary sys- 
tem such as the mass ratio and the unified radius 
ri = Ri/A which normalizes the radius to the 
semi-major axis (actually the accurate ratio can be 
obtained by the LC code). With these values, the 
mean density of the pulsating primary component 
can be precisely determined. From the Kepler's 
law 

4^2^^ 



^"'■•^ G(Mi+M2) 
and the density in solar units 

Mi/i?f 

Mq/R%' 



Pi/Pq 



we get: 



Pi/Pc 



47r2i?3 



MqG{1 + q)rtP^ 



(3) 



(4) 



(5) 



where Porb is the orbit period of the binary system. 

Substituting the values of Porb, 9 and ri from 
the W-D code analysis of the light curve into 
Eq. (5), the mean density is computed, being 
8.47 X 10~^ and 3.67 x 10~^ in solar unit respec- 
tively, which agrees with the density of red giants. 
Then the pulsation constant can be easily calcu- 
lated from: Q — Ppui{pi/ PqY^^ , where Ppui is the 
"pure" pulsation period. The pulsation constant 
Q is an important parameter to tell the intrinsic 
oscillation mode of a pulsating star, especially in 
normal radial mode. The Q value calculated on 
the binary hypothesis, is shown in Table 3 and 
equal to 0.029 and 0.017 for our two candidates, 
respectively. 

The values of rl calculated by the LC code 
in Table 3, 0.47 and 0.44 r espectively, are consis - 
tent with that estimated in lSoszviiski et al.l (|2007l ) 
which notices Ri/A w 0.4 in the binary scenario. 
The result of mode identification shows agreement 
with the theoretica l expectation for red giants 
(|Fox fc WoodI Il982l ) and reveals that 77.7795.29 
is a first overtone pulsating star while 77.8031.42 
is a second overtone pulsating star. The pulsa- 
tion properties derived from a binary model are 
reasonable for a pulsating red giant. 

5. THE EVOLUTIONARY PROPERTY 



A further analysis of the pulsating binary model 
is carried out from the evolutionary property of 
the components, as the stellar properties of the 
red giants are rarely known, by making use of the 
results derived in Section 4. From the 2MASS 
near-infrared photometric data and the empirical 
bolometri c correction factor BC k as a function of 
(J - K) (jBessell fc WoodI Il984[) . and with a dis- 
tance modulus of 18.54 for the LMC, the bolomet- 
ric magnitude can be calculated, and the luminos- 
ity can thus be obtained. Applying the formula 
L = AiraR^T^g, the equilibrium radius of the pri- 
mary star can be computed. Here, Toff is an as- 
sumed value, 3311 K, which is the mean tempera- 
ture for long period variables in LMC. We do not 
calculate the temperature from the color index be- 
cause the secondary component is not understood 
well and its contribution to the observed color in- 
dex is not known. With the mean density and the 
radius, the mass of the primary star is derived. All 
these values are listed in Table 3. 



orb 
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The luminosity of the primary stars in the ta- 
ble is a few thousands of solar luminosity, con- 
sistent with the red giant luminosity. The radii 
of both candidate stars are on the order of hun- 
dreds of solar radius, which are consistent with 
red giants. However, the masses of the two LSPVs 
are too low, and the masses of their components 
would be even lower, according to the mass ra- 
tios in Table 1. Moreover, these low mass stars 
have very high luminosities, which seems impossi- 
ble. For the star 77.7795.29, the mass is 0.31Mq 
and its luminosity is 25O2L0. In the age of uni- 
verse, such a low mass star is impossible to evolve 
up to the red giant branch. For the secondary 
component, the mass is about O.I88M0, while its 
luminosity is 1607Lq, the value of which is com- 
puted from the data of "Li/(Li -I- L2)" in Ta- 
ble 1. It is hard to imagine how such low mass 
stars in a binary system exchange their masses. It 
is difficult to find an evolutionary path to satisfy 
this situation. This problem is also present in the 
star 77.8031.42, with a mass of O.SGM© and a lu- 
minosity of 4839L0, whose component has about 
0.175Mq and 3790^0. If a reasonable mass is ex- 
pected, the density should be upgraded by a factor 
of four or the radius be upgraded by a factor of 1.4. 
Recalling the determination of the density from 
Eq. (5), the density is calculated from the param- 
eters, g, ri and Porb all of which have relatively 
small errors which will not bring about a density 
uncertainty of a factor of four. Besides, the density 
is reasonable for a red giant. On the other hand, 
the radius is derived from the luminosity and effec- 
tive temperature. The luminosity has uncertainty 
from the bolometric correction and the neglect of 
interstellar extinction, and the effective tempera- 
ture also has some uncertainty. If the luminosity is 
higher or the effective temperature lower, the ra- 
dius would be larger and so the mass. However, if 
the mass is indeed higher, the system would then 
have a large separation from Eq.(3) and an even 
more serious problem with the large velocity am- 
plitude. In summary, the pulsating binary model 
requires a low mass incompatible with an evolved 
stellar stage having the observed high luminosity. 

6. DISCUSSION AND CONCLUSION 

The origin of the LSPs is unknown and there 
exist many explanations. The light curves of some 
LSPVs are eclipse-like and it seems that this phe- 



nomenon is due to an invisible component orbiting 
around the pulsating red giant. To test this hy- 
pothesis, we propose a model - pulsating binary, 
and select two LSP stars to analyze their orbit 
motions and pulsation nature. 

On the assumption of binarity, we simulate the 
photometric light curves of the systems by using 
the W-D method. The photometric solutions give 
us a configuration for the binary system: a contact 
system for 77.7795.29 and a semidetached system 
for 77.8031.42 with its Roche lobe fully filled. It 
means that the LSP star may have strong inter- 
action and very probably mass transfer with the 
other component via Roche lobe overflow. How- 
ever, Roche lobe overflow will rule ou t the binary 



hypot hesis if we accept the view of iWood et al 
(|2004[ ). which argues that the mass transfer of 



Roche lobe will result in a short merger timescale 
of about 1000 yrs. Moreover, the calculated ef- 
fective temperatures and the "Li/(Li -I- L2)" val- 
ues in Table 1 suggest that the secondary star 
is also a red giant with a similar temperature 
and luminosity to that of the primary red giant. 
This would lead to a double-lined spectroscopic 
binary. However, no observer of radial veloci- 
ties in these systems has reported seein g spec 
tral li nes from the secondar y star dHinkle et al 



^200i lOlivier fc WoodI 120031: IWood et all 120041: 
iNichoUs et al.l 20091 ). There are also some serious 
problems about the synthesized velocity curves. 
The simulation produces one cycle of velocity 
curve for two cycles of the light curve, while the 
observed velocities show one cycle of the radial ve- 
locity curve for one cycle of the light curve. The 
synthesized full velocity amplitudes for the star 
77.7797.29 and 77.8031.42 are much greater than 
the typical values of other LSPVs. The full ampli- 
tudes of the secondary components are even larger 
and would lead to great velocity separations be- 
tween the primary and the secondary stars, which 
have never been found by spectral observations. 

Fourier analysis is applied to investigate the 
intrinsic oscillation of the LSPVs over the raw 
data sets. Using the parameters {Ri/A and q) 
obtained from the W-D code, the mean densities 
of the LSPVs are deduced and they are consis- 
tent with the red giant phase. Then the pulsa- 
tion constants are obtained by using the classical 
equation Q = Ppuiipi/ PqY^^ to identify the pul- 
sation mode. For star 77.7795.29, one pulsating 
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frequency is detected and it is caused by the first 
overtone radial pulsation, and for star 77.8031.42, 
the only pulsating frequency is caused by the sec- 
ond overtone radial pulsation. These agree with 
the theor etical value for red giants and the con- 
clusion of IWood et all (|l999t ). 

The stellar properties of LSPVs are derived by 
using the information from the pulsating binary 
model. We calculate the bolometric magnitude, 
luminosity, radius and mass of the primary star, 
and find some of them conflict very seriously with 
the evolutionary properties of red giants. In par- 
ticular, the masses for the star 77.7795.29 and 
77.8031.42 are both less than O.4M0, and their lu- 
minosities are too high for such low masses. It is 
difficult to imagine how such low mass stars could 
have such large luminosities. The situation gets 
even worse if we apply the mass ratio and the value 
of + L2)" to calculate the mass and lumi- 

nosity of the secondary star. 

Therefore, the radial velocities and the masses 
computed from the pulsating binary model do not 
agree with some observations and facts about red 
giants. We conclude that the model "pulsating 
binary" has some deficiencies in dealing with the 
observed properties of LSPVs and that the binary 
hypothesis for explaining the LSPs seems unrea- 
sonable. 
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77.8031 .42 













V 















Fig. 1. — Light and velocity curves synthesized by 
the W-D code for 77.7795.29 and 77.8031.42. Also 
shown are the observed R and V band light curves 
from the MACHO project. Red asterisk represents 
the observed data point in R band while the blue 
one is the V band data point. The black solid 
line in the R band panel denotes the theoretical 
synthesized light curve. The bottom panel is the 
synthesized velocity curve. 



Fig. 2. — Power spectrum and light curve of the 
primary star 77.7795.29 in binary system. The 
top five panels show power spectra and the bot- 
tom panel shows the light curve in the R band. 
The black solid light curve is the fit of frequency 
solution described in Sect. 4 of the text. 



S 0.2 r 
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Fig. 3. — The same as Fig. 2, but for the star 
77.8031.42. 
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Table 1 

Photometric solution for stars 77.7795.29 and 77.8031.42. 



77.7795.29 77.80;ii.l2 



classification contact semidetached 

PoA (day) 970.9 1700.9 

Ti (K) 3311(assumed) 3311(assumed) 

T-i (K) 3284 3544±48 

e 0.0016±0.0008 0.0123±0.0164 

i (degree) 62.49±0.27 63.19±2.22 

oj 3.64±0.20 2.91 

q = m2/mi 0.608±0.003 0.486±0.029 

o?i 2.8359±0.0041 2.8748 

(^2 2.8359 2.8943 ±0.055 



Li/(Li +L2) (V band) 0.609 0.561 
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Table 2 

Results of Fourier analysis for stars 77.7795.29 and 77.8031.42. 







77.7795.29 








77.8031.42 






F 


Frequency 


Amplitude 


Phase 


S/N 


Frequency 


Amplitude 


Phase 


S/N 




(c/d) 


(mag) 




(c/d) 


(mag) 




/i 


o.oolio 


0.'2138 


0.3144 


21.6729 


0.0011 


0.1464 


0.46905 


19.3482 


/2 


0.0100 


0.0529 


0.6576 


22.4815 


0.0112 


0.0770 


0.7190 


44.7148 


/3 


0.0162 


0.0327 


0.0015 


13.9081 


0.0022 


0.0613 


0.2773 


35.6114 


/4 


0.0011 


0.0363 


0.9538 


15.4348 


0.0116 


0.0446 


0.7607 


25.9384 


h 


0.0041 


0.0280 


0.6792 


11.9004 


0.0120 


0.0404 


0.6926 


23.4930 


h 


0.0130 


0.0205 


0.8848 


8.7069 


0.0040 


0.0260 


0.1102 


15.1249 


fr 


0.0174 


0.0189 


0.4320 


8.1379 


0.0162 


0.0267 


0.1793 


15.4983 


fs 


0.0031 


0.0199 


0.0318 


8.0322 


0.0029 


0.0351 


0.8120 


20.3622 




0.0049 


0.0199 


0.6315 


8.4640 


0.0169 


0.0224 


0.8329 


13.0511 


JlO 










0.0108 


0.0311 


0.7704 


18.1114 












0.0070 


0.0287 


0.0900 


16.6748 


J12 










0.0057 


0.0236 


0.3752 


13.7059 


/l3 










0.0065 


0.0242 


0.2452 


14.0897 


/14 










0.0006 


0.0206 


0.8211 


11.9582 
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Table 3 

Parameters of the primary stars 



Star 


Frequency 
(c/d) 


Ri/A 




Q 


Mode 


Mboi 
(mag) 


L/Lq 


R/Rq 


M/Mq 


77.7795.29 


0.0100 


0.47 


8.47E-8 


0.029 


IH 


-3.75 


2502 


153 


0.31 


77.8031.42 


0.0112 


0.44 


3.66E-8 


0.017 


2H 


-4.46 


4839 


213 


0.36 
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